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We investigate the effect of drag force on the enstrophy cascade of two-dimensional Navier-Stokes 
turbulence. We find a power law decrease of the energy wave number (k) spectrum that is faster than the 
classical (no-drag) prediction of k~?. It is shown that the enstrophy cascade with drag can be analyzed 
by making use of a previous theory for finite lifetime passive scalars advected by a Lagrangian chaotic 
fluid flow. Using this we relate the power law exponent of the energy wave number spectrum to the 
distribution of finite time Lyapunov exponents and the drag coefficient. 


PACS numbers: 47.27.Eq, 47.52.+j, 83.10.Ji, 83.50.Ws 


Two-dimensional Navier-Stokes fluid turbulence has 
been of interest for many years. Examples where two- 
dimensional (2D) Navier-Stokes turbulence is potentially 
relevant include soap film flow [1-6], rotating fluids [7], 
magnetically forced stratified fluids [8,9], plasma in the 
equatorial ionosphere [10], and the Earth’s large scale 
(>500 km) atmospheric dynamics [11]. In these situations 
there are regimes where drag is an important physical 
effect [12]. In this case the describing Navier-Stokes 
momentum equation is 
= tv: VW = —p !Vp + vı V?v — vov +f, (1) 
where vo is the drag coefficient, vı is the viscosity, f is an 
external forcing term, and incompressibility (V - v = 0) 
will be assumed. If the forcing is localized in k space 
with a characteristic wave number of kp, it is well known 
that, in the absence of drag (vp = 0), an energy cascade 
to larger scales, E(k) ~ k~>/3 (k < ky), and an enstro- 
phy cascade to smaller scales, E(k) ~ k~* (k > Kp), are 
expected [13,14]. Our concern in this paper will be the 
effect of drag (the term — vov) on the enstrophy cascade 
in two-dimensional turbulence forced at low wave num- 
ber (long wavelength). In the absence of drag (vo = 0), 
many experimental and numerical results are consistent 
with the theoretical prediction of a k~* dependence of the 
energy wave number spectrum E(k) at wavelengths below 
the forcing wavelength [5,9,15] [while some other experi- 
mental and numerical results show that E(k) ~ k~> be- 
havior is not conclusive yet [1—4,16] ]. In the presence of 
drag, we find that there is still a power law dependence of 
the energy spectrum on k, but the power law exponent is 
greater than three, 


E(k) ~ k~8t8), ¿>0. 


This might be relevant to results of experiments where 
energy spectra steeper than k7? were reported [1-4,16]. 
Furthermore, we find that the exponent increase é can be 
quantitatively related to the drag coefficient vo and the 
chaotic straining induced by the flow. The key point in 
obtaining a theory for the exponent increase é is our argu- 
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ment (given subsequently) that the presence of drag funda- 
mentally changes the character of the enstrophy cascade, 
making it possible to analyze the high wave number com- 
ponents by the use of a recent theory for the wave number 
spectrum of a finite lifetime passive scalar in a Lagrangian 
chaotic fluid flow. In particular, we will argue that the high 
k components of vorticity behave like a passive scalar ad- 
vected by the Lagrangian chaotic low k components of the 
velocity field. 

In our studies, we numerically solved the vorticity equa- 
tion obtained by taking the curl of Eq. (1); 

ðw 


A A 


z : Vo = (-1)"*!y,V"@ — vow + So(x), 


(2) 
withy = 2 X Vw, Vy = w, where y is the stream func- 
tion, œ = 2+ V X vis the vorticity, and we have replaced 
the term v,V*w@ resulting from (1) by a hyperviscous 
damping, (—1)"*!v,V?"@. Se is a source function of 
the vorticity (= Z- V X f). Hyperviscosity (n = 4 for 
our numerical computations) is a commonly used device 
to enhance the wave number power law scaling range in 
numerical turbulence simulations. Our simulations employ 
a two-dimensional square domain, [—7, m] X [—7, 7], 
with periodic boundary conditions, and S(x,y) = sin2y 
is used for the source function of the vorticity. All the 
simulation results presented here use a spatial grid of 
1024? and v4 = 5 X 107”. 

The energy spectrum resulting from a simulation with 
vo = 0.1 is shown in Fig. 1. The energy spectrum is ob- 
tained by time averaging of instantaneous energy spectra 
measured at every 5 time units over a time duration of 40. 
The dashed line corresponds to E(k) ~ k~>°, where the 
exponent is from our theory which will be explained later. 
A clear power law behavior over almost one decade in k 
space is observable until the spectrum rolls off at high k’s 
due to the hyperviscosity. Figure 2 shows plots of k* E(k) 
versus k for three cases with vp = 0,0.1, and 0.2 applied 
at high k (k > 6), but with the same drag, vp = 0.1, ap- 
plied at low k (k = 6) in all three cases. The upper- 
most curve (circles) is for a high k drag of vo = 0, the 
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FIG. 1. Energy spectrum with a drag force coefficient vo = 
0.1. The slope of the dashed line is —3.5 from the theory based 
on G(h). 


middle curve is for vọ = 0.1, and the lowermost curve is 
for vo = 0.2. Applying the same low k drag allows us to 
compare effects of drag for high k’s while keeping similar 
large scale dynamics of the flows [17]. This is impor- 
tant since we will show that, with drag, the major contri- 
bution to the straining of high k vorticity components is 
from the longest length scale (lowest k) flow components. 
While the energy spectrum for our vọ = 0.0 case shows 
~k~> behavior, consistent with the classical Batchelor- 
Kraichnan prediction, the energy spectrum shows steeper 
power law behavior as we increase the drag coefficient. A 
summary of the simulation results appears in Table I. 

In order to start to formulate a theory of E(k) in the 
presence of drag, we first consider the straining rate, n (r), 
on eddies of size r, 


= / /r 
Pee (ive +r) rer a ia ElK? dk. 


r? 7 

(3) 
where, in setting the lower limit of the integration over 
k at kp, we have in mind the situation of our numerical 
simulation where ky ' is of the order of the system size 
and there is appropriate long wavelength damping [17]. In 
(3) the contribution of E(k) from high wave numbers be- 
comes negligible for large k (small r) if E(k) ~ k~S*8) 
with € > 0 (which is the case when drag is present). For 
the classical case without a drag term where é = 0, the 
above quantity diverges logarithmically at large k (small r) 
and this leads to Kraichnan’s prediction of a weak logarith- 
mic correction to the power law energy spectrum, E(k) ~ 
k~? log™!/ (kL) [18]. In contrast, if E(k) ~ k~S+® with 
é > 0, then for small r the straining rate becomes con- 
stant, ie " E(k)k? dk = fo E(k)k? dk, indicating that the 
straining of smaller scales comes from the larger scales. 
Thus, for the € > 0 case, it is a good approximation to 
assume that the vorticity components at sufficiently high 
k are passively advected by the chaotic velocity field of 
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FIG. 2. Energy spectrum with various drag coefficients vo. 
The circles, the crosses, and the stars are for vo = 0.0, vp = 0.1, 
and vp = 0.2, respectively. 


large structures and this approximation becomes better and 
better as k is increased. In treating the case with drag 
we shall use Lyapunov exponents. We thus need to dis- 
cuss the relevance of Lyapunov exponents to our prob- 
lem. Crudely, we estimate the Lyapunov exponent h as the 
straining rate at vanishingly small separation, i.e., r = 0 
in Eq. (3), h ~ Lhe, k? E(k) dk]'/. In the presence of vis- 
cosity (or hyperviscosity), the resulting high k cutoff to 
E(k) implies a finite Lyapunov exponent. However, de- 
pending on the exponent in the power law dependence of 
E(k), different wave numbers, i.e., different scale lengths 
will determine the value of the Lyapunov exponent. With 
no drag we expect E(k) ~ k~? and the Lyapunov exponent 
is dependent on contributions from flow components with 
wave numbers down to the viscous cutoff scale. [This is 
signaled by the logarithmic dependence of f ie " E(k)k? dk 
on r.] This means that the Lyapunov exponents only de- 
scribe straining of scales smaller than the viscous cutoff 
scale, and hence are not relevant to the enstrophy cascade. 
On the other hand, if E(k) ~ k~S*® with é > 0, as we 
claim for the case with drag, then the situation is different. 
When é > 0, Lyapunov exponents are relevant in the ab- 
sence of a high k cutoff and are determined by the largest 
scale flow components. Thus, for € > 0, the Lyapunov 
exponents describe the straining of scales small compared 
to kp t which, however, may still be large compared to a 
viscous cutoff. 


TABLE I. Comparison of the power law exponents from the 
direct numerical simulations, (3 + éns), to the theoretical re- 
sults, (3 + én). The spectral ranges linear fitted to measure the 
power law exponents are log,)k = [1.0, 2.0] for the direct nu- 
merical simulation. 











volk > 6) 3 + Ens 3 + Én 
0.0 3.1 3.0 
0.1 3.6 3.5 
0.2 4.1 4.0 
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The discussion in the previous paragraph implies that the 
situation for w at high k is the same as the situation for a 
passive scalar œ convected by a Lagrangian chaotic flow v 
in the presence of a loss term such that ¢ obeys an equa- 
tion [analogous to (2)] of the form 0¢/dt + v : Vd = 
DV? ¢ — vod + S(x), where —voġ is the loss term. 
(This type of spectral passive scalar problem is of interest 
for certain reacting chemical scalars and also for plankton 
concentration in the ocean [19—21].) The wave number 
spectrum for such a passive scalar has recently been con- 
sidered in Ref. [20], and we can make direct use of that 
reference. According to Ref. [20], the wave number power 
spectrum of ¢ behaves as Ey(k) ~ kCt) where é is de- 
termined by the probability distribution function P (h, t) for 
finite time (t) Lyapunov exponents (h). In order to define 
this function, consider a disk of fluid of infinitesimal ra- 
dius 67ro originally located at a point r, and evolve the disk 
forward for a time t. At time t the disk has been deformed 
into a differential ellipse whose major radius is ôr,. The 
time ¢ finite time Lyapunov exponent for the initial point 
r is h(r,t) = In(6r,/dro). P(h,t)dh is the probability 
that A(r, t) is between h and h + dh given that the initial 
point r is chosen randomly. For large t, it can be shown 
that InP (h, t) has an asymptotic form InP (h, t) ~ —tG(h), 
where G (h) is concave up [G”(h) > 0] and has a minimum 
at h [i.e., G'(h) = 0] with G(h) = 0, where A is the usual 
infinite time Lyapunov exponent [22]. According to [20], 
the spectral exponent is determined from G (A), 


& = mini[G(h)/h] + (2vo/h)}. (4) 


As argued above, at k >> ky we can treat w as a passive 
scalar in a Lagrangian chaotic fluid flow so that E, ~ 
k-"+8) Also, since V: v =0, w =2: V Xv, we 
have E(k) = k?E,(k), and E(k) ~ k-S*® with é 
given by (4). 

We test the applicability of this theory to the enstrophy 
cascade with drag by numerically obtaining G(h) from a 
histogram approximation to the finite-time Lyapunov ex- 
ponent distribution P(h, t). The latter is obtained by com- 
puting h for 4 X 104 initial conditions spread uniformly 
in space. A cubic polynomial fit to +7! InP (h, t) is utilized 
to obtain G(h). Theoretical power law exponents obtained 
from Eq. (4) and the numerical G(h) are compared with the 
results from full numerical simulations in Table I. Within 
the accuracy of our computations reasonable agreement is 
obtained. 

It is also of interest to consider the possible power law 
scaling of the vorticity structure functions, 


(la(x +r) — olx) ~ r“, (5) 


for small r = |r| and r > (viscous cutoff length). Re- 
cently the structure functions for Lagrangian chaotic 
advection of finite lifetime passive scalars have been 
considered [23]. By our previous arguments we can 
directly apply the result of Ref. [23] to the vorticity 
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structure functions for two-dimensional turbulence with 
drag, 


éq = mimg,[G(h)/h] + (qvo/h)}, (6) 


which generalizes (4). That is, 2) agrees with £ in (4). 

For the case without drag (vo = 0), Eq. (6) indicates 
that €, = 0 for all q [recall that G(h) = 0]. This indicates 
the absence of small scale intermittency in the enstrophy 
cascade for vp = 0, in agreement with previous work (see 
[9] and references therein). This situation is in striking 
contrast with three-dimensional turbulence where intermit- 
tency is clearly present. Equation (6), however, shows the 
important result that drag (vo > 0) leads to intermittency 
also in the two-dimensional case (¢, # 0 and not propor- 
tional to q). 

In conclusion, we have considered the effect of linear 
drag on the enstrophy cascade of two-dimensional Navier- 
Stokes turbulence. We find that, as in the case without 
drag, enstrophy cascade results in power law behavior of 
the energy spectrum. However, the power law exponent is 
different from the classical Batchelor-Kraichnan prediction 
(—3) and is determined by the distribution of finite time 
Lyapunov exponents and the drag coefficient. 
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